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Abstract 

For systems that involve particle production through branching pro- 
cesses the concept of chaos is explored. The measures that can describe 
their behaviors are investigated. Monte Carlo simulation is used to gen- 
erate events according to perturbative QCD and an Abelian model. It is 
shown how the measures proposed distinguish the two cases in ways that 
characterize the chaotic behavior. 



It has been known for some time that the nonlinear, non- Abelian dynamics of the 
classical Yang-Mills held has chaotic solutions [1, 2]. More recently, it has been shown 
by lattice calculation that the classical non- Abelian gauge theory generally exhibits 
deterministic chaos and that the Lyapunov exponents can be numerically determined 
[3, 4, 5]. How to extend the investigation to the quantum theory is, however, unclear 
inasmuch as the notion of quantum chaos for such dynamics is not well defined [6, 7]. 
In this paper we take the hrst step in that direction, not just going into the quantum 
dynamics of a nonlinear theory, but into the realm of particle production of quantized 
fields. 

In this uncharted territory we have very little guidance on what to study in search 
for signs of chaos. It is not clear what a trajectory is in QCD, even less the distance 
between two trajectories. What exactly is the Kolmogorov entropy, well defined in 
classical dynamics [8], is also unclear in the multiparticle production problem. Our 
hrst objective is therefore to the find some measure that can play the role of distance 
between trajectories and some other quantity that conveys the loss of information in 
the final state. 
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In order to know whether or not we have found the right measures, it is necessary 
to test our ideas on different dynamical problems and show that the measures can 
effectively distinguish different characteristics. We shall therefore concentrate on 
two branching processes in particle production. One is the pure gluon theory in 
perturbative QCD; the other is a cascade model (to be referred to as the x model) 
that has none of the features of the gauge theory. We shall regard the gluons (partons) 
as particles and ignore hadronization so that we may focus on the issue of chaoticity 
in the branching dynamics. 

In either case we start with a parton having high time-like virtuality Q 2 and use 
computer simulation to study the evolution of the system through branching. The 
Monte Carlo code generates the parton momentum distributions that in the case of 
QCD satisfy the Altarelli-Parisi evolution equations [9, 10, 11]. The splitting functions 
for the two cases are drastically different, and the ways the virtualities degrade are 
also different. They give rise to the diverse behaviors that emerge. The question is 
how to find a quantitative measure of the diversity that is useful from the point of 
view of chaoticity. 

One of the difficulties in the problem posed is that time is not an explicit variable 
when calculating momentum distribution. Moreover, with the number of degrees of 
freedom increasing in the cascade process, there is no obvious way to generalize the 
definition of a trajectory in classical dynamics. However, branching has an irrefutable 
physical notion of the direction of time, and parton multiplicity n increases with time 
so long as only tree diagrams (without recombination) are considered, which we do. 
Thus we take two conjugate views: the "temporal" development (parameterized by n) 
without regard to the parton momenta, and the "geometrical" characteristics in the 
momentum space at the end of the evolution. In some rough sense this corresponds 
to the two views of the classical systems that are related by ergodicity. 

Concerning trajectory and distance between trajectories, our thinking goes as 
follows. The classical field intensity is replaced in the quantum case by the number 
n of quanta, and the distance between two field configurations can be represented 
by the variance D 2 of the n distribution, P n . This distribution arises after repeated 
simulation from the same initial condition, i.e., a fixed initial virtuality Q 2 . Quantum 
fluctuation is enough to replace the small variation of the initial condition in classical 
dynamics. With the average (n) regarded as a measure of the evolution time, D 2 as a 
function of (n) can then be regarded as the quantum analog of the classical distance 
function dit). This (n) is the average multiplicity, not just in the final state, but 
along the evolution process. 

More precisely, let us focus on one event and use a tree to represent a particular 
branching process. Regardless of the virtuality q 2 of any line, all vertices of the same 
generation are put at the same branching level and a label i is given to each generation 
between two levels, starting with i = for the initial parton at Q 2 . By ignoring q 2 
in the genealogy of the tree, we are emphasizing topology and overlooking momenta. 
Let rii denote the number of partons at the ith generation. We shall use the event 
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averaged in lieu of time, although a simple linear or exponential relationship 
between the two is not expected. In fact, may not increase monotonically with 
z, even if n 8 - does; in such a case we use only the increasing portion. For the measure 
that plays the role of distance, we define the normalized variance 



Vi= {{nl) ~ (n t ) 2 ) / (n t ) 2 (1) 

This differs from the second cumulant moment A 2 by a term — I/(n 8 ), and is better 
than it because A 2 can become negative (thus precluding log-log plots), while V{ is 
always positive. A rapid increase of V{ with can be interpreted as the analog of 
the "divergence of nearby trajectories" in classical dynamics. Clearly, if V{ remains 
constant or decreases, one does not become more ignorant about the state of the 
system as it evolves, and thus it cannot be regarded as exhibiting chaotic behavior. 

For MC simulation we follow the procedure described by Odorico [10]. For both 
pure-gauge QCD and the \ model we start with virtuality Q 2 and end at Qq. The 
splitting function for g ^ gg is 



P(z) = 6 



1 — z z 
z I — z 



(for QCD) 



where z is the momentum fraction of the daughter parton in the frame where the 
mother parton's momentum is I. In terms of the Sudakov form factor S(Q 2 } Qq) 
one can calculate the probability of emitting a resolvable gluon occurring between 
Q 2 and Q^; when it occurs in the simulation at q 2 } then z is chosen to be in the 
range z < z < I — z , where z = Ql/q 2 - The daugthers evolve separately from 
the maximum virtualities, q\ = z q 2 and q\ = (I — z)q 2 } according to the same 
procedure as applied to the mother. The running coupling constant is, as usual, 
a s (q 2 ) = 47r/[II log(g 2 /A 2 )], where we have set N c = 3, Nj = 0, and shall use 
A = 250 MeV. Branching ceases when q 2 reaches or below. 
In the x model we use the splitting function 

P(z) = 6z (1 — z) (for the x Model) (3) 

There is no infrared or collinear divergence and therefore no evolution. Nevertheless, 
we introduce Q 2 dependence by requiring that the daughter virtualities be zq 2 and 
(I — z)q 2 } when the mother virtuality is q 2 } and z can be any value between and 1. 
We require branching to occur successively until the virtualities of all lines reach Qq. 
We consider this model because it examplifies the Abelian dynamics without infrared 
and collinear divergences. Yet the multiplicity of particles produced depends on Q 2 
so that the result of branching can be compared with the non- Abelian case, although 
for very different Q/Qo. 

We have simulated these two branching processes by running I0 5 events each, 
using Qo = 1 GeV. The results on V{ vs are shown in the log-log plot in Figure I 
for various values of Q /Qo indicated. While it is hard to produce high multiplicity in 
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QCD unless Q/Qo is extremely large, particles are copiously produced in the \ model 
at moderate Q/Qo. The general features of Vi vs for the two cases are markedly 
different. The QCD result shows a power-law increase in the high range 



V oc (n t ) 



0.4 



(4) 



where the exponent n is approximately independent of Q 2 . The \ model, on the other 
hand, shows a rapid rise initially, but followed by a precipitous drop after reaching 
a maximum. Indeed, the maximum Vi decreases with increasing Q 2 . Clearly, this is 
not that suggest chaotic behavior. 

For gluon branching the monotonic increase of Vi with implies that, as the 
branching proceeds, how many particles are produced in any event becomes more and 
more unpredictable. The power-law dependence may be regarded as the analog of the 
exponential increase with time of the distance between classical trajectories that are 
initially close by. However, there is no way to relate n to the Lyapunov exponents, 
since among other differences the notion of time is not well defined here. By itself 
k ~ 0.4 does not indicate how chaotic the behavior is. There is a need for another 
measure of chaoticity. 

If Figure I is viewed as the analog of the description of the temporal behavior, 
another place to search for signs of chaos is in the phase space of the particles. As 
the system evolves, more and more information is lost on where the partons are (or 
more precisely, what their momenta are), so entropy increases not only because of the 
increase of the number of particles, but also because of the dynamical fluctuations in 
their momenta. We therefore consider a multifractal description of that fluctuation 
and focus on the information dimension as a characterization of the entropy of the 
system [8]. 

Since at each vertex of branching a daughter's momentum fraction z is known 
in the simulation, the momentum x of a final particle as a fraction of the initial 
particle is therefore calculable. It is x = Yli where Zi is the momentum fraction of 
the descendant at the ith generation. Since the particle distribution p(x) is highly 
peaked near x = 0, it is smoother to examine the distribution in the cumulative 
variable X, defined by [12] 



where x\ and x 2 are two extreme points in the distribution p(x), between which X 
varies from to I. In terms of X the inclusive distribution p(X) is constant. For each 
event the fluctuation in X space is then studied by dividing the interval < X < I 
into M bins and calculating the factorial moments 




(5) 



M 



f q (M) = M- 1 £ n 3 {n 3 - I) • • • {n 3 - q + I) 



(6) 
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where nj is the multiplicity in the jth bin. After averaging over all events, the 
normalized factorial moment 

F q = (f q )/(fi) q (7) 

is known to contain no statistical fluctuations [13]. 

What interests us is its behavior near q = 1, where we can extract the information 
dimension D\. To that end it is necessary to extend the definition of F q in (6) to 
noninteger q. A method for achieving that without losing the attribute F q = 1 for 
Poissonian fluctuation has recently been developed [14]. Using that method we have 
calculated F q for a continuous range of q, as shown in Figure 2, for both QCD and 
the x model. The results for the two cases are very different and provide a distinct 
contrast between them. They are in accord with the temporal behavior shown in 
Figure 1 in that, for q > 1, F q < 1 in the \ model, meaning that the distribution is 
sub-Poissonian, while F q > 1 in QCD, indicating large fluctuations. The geometrical 
properties in the X space are not revealed until we investigate the M dependence. We 
find that in both cases F q are not sensitive to M , as can be seen from the various lines 
in Figure 2 corresponding to different M values. Thus the behavior has no interesting 
multifractal property: D\ ~ 1 in both cases. 

The origin of the lack of significant M dependence can be traced to F q itself, 
where the event averaging cancels out the fluctuations. Event by event the values 
of Fy = f q /(fi) q } where e labels an event, fluctuate greatly, especially when M is 
large. To quantify the degree of that fluctuation we define event-averaged (vertical) 
moments of the (horizontal) _F e moments 

C p , q (M) = (F>(M)) / (F q (M))> , (8) 

where (F^ = A/" -1 J2 e (Fq Y \ N being the total number of events. We then calculate 
C Pj q(M) for < p < 2, and q = 2,3,4. It is found that C Pjq (M) indeed exhibits 
significant M dependences, as shown in Figure 3. In all cases of p and q, the \ 
model has smaller C Pjq compared to QCD, implying smaller fluctuations of _F e . The 
M dependences do not show linearity over any extended range in the log-log plots, 
the best being from M = 5 to 20. In that range we write 

C P JM) oc M^W . (9) 

From the slope ip q (p) in the neighborhood of p = 1, we can calculate the index /i g , 
defined by 

N = ^q(p) \p =1 ■ ( 10 ) 

The result of our simulation yields the values fi q = 0.0061, 0.054, 0.23 for q = 2, 3, 4 in 
the case of QCD, and 0.0014, 0.010, 0.046 in the x model, respectively. Clearly, for 
each 5, //(Q CD ) is significantly larger than . We now give a physical interpretation 
of fi q as an entropy index to be used as a new measure of the event fluctuation in 
branching processes. 



5 



If we define P g e = F^j Y, e Fqi an d then define H Pjq = J2e(Pq) P 7 we have an entropy 
in the event space 

S, = -£P;inP; = -^lnff p | p=1 . (11) 
H Pj q and C Pjq can be related by their definitions, yielding 

-^-lnC p ? \ p=1 = InjV - S q . (12) 
dp 

From (9) and (10) we obtain 

Sq=\n{MM-^) , (13) 

apart from a possible additive term independent of Af and M. fi q appears to be 
related to the information dimension, but it is not because S q is not the entropy 
defined in the X space, which is the one that is divided into M cells. The event 
space in which S q is defined has not been partitioned into small cells. The meaning 
of (13) can be seen in two extreme cases: (a) if F^ is the same for every event, then 
Pq = 1/A/", and S q = bxV ; (b) if only one event has _F e ^ 0, and _F e = in all others, 
then S q = 0. Thus case (b) is more ordered in the event space than (a); that is, it is 
more disordered to spread out an observable (F^ here) over all events (even if F^ = 
constant) than to confine it to a few events having nonzero values (analogous to the 
increase of entropy of an expanding gas). [The case of all F^ = is excluded from 
consideration in order to render P g e meaningful.] Thus S q decreases when there are 
more events with _F e = 0, signifying more order in the event space. From (13) we see 
that fiq is a measure that decrease, which in turn implies more fluctuation in F^. 

At large M only large spikes in small bins contribute to F^ } especially when q is 
large. Events with large spikes are rare. Consequently, the fluctuation in F^ from 
event to event becomes more pronounced with increasing q. That behavior is now 
quantified by fi q . We may therefore use fi q to characterize the spatial properties 
of the chaotic behavior of a branching process. We have, however, at this stage no 
quantitative criterion on how small fi q must be in order to signify no chaotic behavior. 

We can relate the classical and quantum problems in our description in the event 
space as follows. Consider a neighborhood N e around an initial point in phase space 
for a classical trajectory. For chaotic dynamics, starting the system from any point 
in N e leads to widely different trajectories at sufficiently long time later. We may 
regard Af of these trajectories in N t as corresponding to AT branching events all with 
the same initial virtuality, but having different outcomes. F^(M) describes the final 
state of the system for the eth event, and fi q describes the degree of fluctuation of F^ 
from event to event. Sufficiently large values for the index fi q therefore signify chaotic 
behavior of the branching process. 

In conclusion, we have found features about QCD branching that are not shared 
by the \ model, which represents Abelian branching. Because of the non-classical 



6 



nature of the system, we have had to search for new measures and observables. The 
dependence of V{ on reveals the "temporal" behavior, while F q , C Pjg , and fi q 
describe the "spatial" properties. All these measures taken together give a collective 
description of the degree of chaoticity in a branching process. We have found that V{ 
increases with in QCD, while it decreases for the \ model. The dependences of F q 
on q are totally different for the two much larger for QCD than 

for the x model. These results collectively suggest that QCD branching is chaotic, 
while the \ model is not. Among the measures considered, V{ (and higher moments of 
rii, which could also have been considered) vs contain detailed information about 
the branching process from generation to generation, but they are not accessible to 
the experiment. The others describe the characteristics of the final state, and can be 
determined experimentally in most high-energy collisions. The entropy index fi q is 
most unusual and deserves further investigation. 

One of us (RCH) is grateful to K. Geiger, B. L. Hao, S. G. Matinyan, and B. Miiller 
for helpful discussions. This work was supported in part by the U.S. Department of 
Energy under Grant No. DE-FG06-9f ER40637. 
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Figure 1: Normalized variance vs average multiplicity at various generations in the 
branching processes. 
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Figure 2: Factorial moments of continuous order for various bin sizes. Q/Qo = 10 3 
for QCD and 20 for the x-model. 
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Figure 3: C Pjq vs M for various values of p and 5. Q/Qo = 10 3 for QCD and 20 for 
the x-model. 
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